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Abstract 

Ecological factors often shape demography through multiple mechanisms, mak- 
ing it difficult to identify the sources of demographic variation. In particular, 
conspecific density can influence both the strength of competition and the pre- 
dation rate, but density-dependent competition has received more attention, 
particularly among terrestrial vertebrates and in island populations. A better 
understanding of how both competition and predation contribute to density- 
dependent variation in fecundity can be gained by partitioning the effects of 
density on offspring number from its effects on reproductive failure, while also 
evaluating how biotic and abiotic factors jointly shape demography. We exam- 
ined the effects of population density and precipitation on fecundity, nest sur- 
vival, and adult survival in an insular population of orange-crowned warblers 
{Oreothlypis celata) that breeds at high densities and exhibits a suite of traits 
suggesting strong intraspecific competition. Breeding density had a negative 
influence on fecundity, but it acted by increasing the probability of reproduc- 
tive failure through nest predation, rather than through competition, which was 
predicted to reduce the number of offspring produced by successful individuals. 
Our results demonstrate that density-dependent nest predation can underlie the 
relationship between population density and fecundity even in a high-density, 
insular population where intraspecific competition should be strong. 
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Introduction 

Quantifying the demographic effects of population density 
is essential for understanding how populations are regu- 
lated (Murdoch 1994), but poses a major research chal- 
lenge because population density can act via multiple 
mechanisms and in concert with other ecological factors 
(Ostfeld and Canham 1995; Hixon and Carr 1997; Forch- 
hammer et al. 1998). For example, population density can 



shape fecundity through increased resource competition 
(Rodenhouse et al. 2003) or increased density-dependent 
offspring mortality (Arcese et al. 1992). Density-depen- 
dent reductions in survival can act on both juveniles 
(Clutton-Brock et al. 1987; Gaillard et al. 1998; Harms 
et al. 2000) and adults (Forrester 1995; Frederiksen and 
Bregnballe 2000) through altered rates of starvation, pre- 
dation, parasitism, or disease (Dempster 1971; Walde and 
Murdoch 1988; Hochachka and Dhondt 2000; Holbrook 
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and Schmitt 2002). Other ecological factors, such as food 
abundance and weather, exhibit the same potential to affect 
demography via multiple mechanisms, and also mediate 
the effects of population density on vital rates (Coulson 
et al. 2001; Sillett et al. 2004). Therefore, developing a 
mechanistic understanding of demographic variation 
requires partitioning the multiple processes by which eco- 
logical factors shape patterns of survival and reproduction. 

The mechanisms by which population density influ- 
ences fecundity are often obscured in natural populations, 
in part because fecundity is comprised of several compo- 
nents. Fecundity is defined as the number of young pro- 
duced over a breeding season, and variation in this 
demographic rate arises because (1) many individuals fail 
to reproduce successfully and (2) successful individuals 
vary in the number of young they produce (Clutton- 
Brock 1988; Newton 1998). However, we rarely know 
whether the ecological processes that cause reproductive 
failure, here defined as producing no young over a breed- 
ing season, are distinct from those affecting offspring 
number. For example, studies of birds show that while 
food availability is often the most important factor limit- 
ing the number of offspring produced by successful indi- 
viduals (Lack 1947; Martin 1987; Godfray et al. 1991; 
Nagy and Holmes 2005; Sofaer et al. 2013), low food 
abundance can also lead to failure by causing individuals 
to forgo breeding (Southern 1970; Grant et al. 2000; 
Jenouvrier et al. 2003; Langin et al. 2009). Similarly, nest 
predation is the primary cause of avian nest and repro- 
ductive failure (Ricklefs 1969), but can also affect off- 
spring number by influencing clutch size (Skutch 1949; 
Martin et al. 2000; Eggers et al. 2006) or the number of 
young fledged (Zanette et al. 2011). Because population 
density can alter food availability (Arcese and Smith 1988; 
Newton 1998) and nest predation risk (Gunnarsson and 
Elmberg 2008), it can act via either factor to affect off- 
spring number or the probability of reproductive failure. 
Identifying and partitioning the processes that affect each 
component of fecundity can therefore reveal the mecha- 
nism through which density dependence shapes demo- 
graphic variation. 

Island populations have been used as model systems to 
investigate the mechanistic basis of demographic varia- 
tion. These studies typically highlight the importance of 
both competition (Blondel 2000; Brouwer et al. 2009; 
Nevoux et al. 2011) and the abiotic factors affecting food 
abundance (Grant et al. 2000). Intraspecific competition 
is thought to be strong in island populations because they 
occur at high densities relative to their mainland counter- 
parts and are often characterized as having fewer preda- 
tors (MacArthur et al. 1972; Yeaton 1974; Adler and 
Levins 1994). However, despite the past emphasis on den- 
sity-dependent competition, nest survival (the probability 
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a nest fledges one or more young) in island bird 
populations can also be density dependent (Arcese et al. 
1992), and nest predation rates on islands may be no 
lower than those on mainlands (Covas 2012). These pat- 
terns emphasize the need to determine how different 
mechanisms of density dependence and the abiotic drivers 
of food abundance jointly shape the demography of 
insular populations. 

Here, we analyze the ecological correlates of demo- 
graphic variation in an island population of orange- 
crowned warblers [Oreothlypis celata; Fig. 1). Our study 
population breeds at a higher density than mainland pop- 
ulations (see below) and exhibits traits associated with 
increased intraspecific competition (e.g., high testosterone, 
Horton et al. 2010; high aggression, Yoon et al. 2012). We 
therefore expected density-dependent competition to 
underlie demographic variation in this system, along with 
variation in rainfall, which is a driver of food abundance 
(Morrison and Bolger 2002; Sofaer et al. 2013). We stud- 
ied the effects of population density and rainfall on two 
major vital rates: fecundity and annual adult survival. Our 
analysis of fecundity separated the processes affecting the 
offspring number of successful individuals from those 
affecting the probability of reproductive failure, allowing 
us to assess the demographic effects of intraspecific com- 
petition and nest predation. By partitioning the effects of 
different processes on a single demographic rate, we were 
able to clarify the mechanisms of density dependence and 
challenge our assumptions about the relative importance 
of competition and predation. 

Methods 



Field methods 

We studied breeding orange-crowned warblers in Bulrush 
Canyon (33°20'N, 118°26'W) on Santa Catalina Island, 




Figure 1. Orange-crowned warbler {Oreothlypis celata) on Santa 
Catalina Island, California. Photo by Dr. Moo-Boo Yoon. 
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CA (hereafter, Catalina Island), each spring from 2003 to 
2009. Our study population belongs to the O. c. sordida 
subspecies, which is endemic to the California Channel 
Islands and remnant patches of coastal chaparral habitat 
on the adjacent mainland. Most adults within a 7-ha 
study area were captured in mist-nets using conspecific 
playback and given a unique combination of colored plas- 
tic leg bands and a numbered US Geological Survey alu- 
minum band. This study area represented the central core 
of our study plot, which was intensively surveyed in all 
years. Each March, surveys were conducted throughout 
the study area and adjoining habitat to resight color- 
banded individuals for analyses of survival. We observed 
breeding pairs every 1-3 days throughout each breeding 
season to map territory boundaries relative to a 
25 X 25 m grid system, and to locate and monitor all 
nesting attempts, including renesting and double brood- 
ing (see Sofaer et al. 2013). The number of territorial 
pairs in our study area ranged from 23 (3.4 territories/ha) 
to 43 (6.3 territories/ha), and for each pair, we calculated 
the distance to the nearest neighbor (m) based on the 
center of each territory polygon, as defined in ArcMap 
(ESRI, Inc., Redlands, CA). 

Fecundity and nest survival 

We quantified the annual fecundity of 181 territorial war- 
bler pairs. In order to exclude pairs that may have had an 
undetected renesting attempt, we limited our fecundity 
analysis to pairs that either successfully fledged young at 
any time during the breeding season or had an active, 
monitored nest after April 15th; birds that failed after this 
date were unlikely to successfully renest. However, our 
analysis of fecundity did include closely monitored pairs 
without any active nests in 2007, when an extreme 
drought caused nearly all pairs to skip breeding and only 
one monitored pair laid eggs (Langin et al. 2009); no data 
from that year were included in our analysis of nest sur- 
vival. The average number of pairs included in our fecun- 
dity analysis each year was 26 (range: 13-33 pairs; Table 
SI). We calculated annual fecundity of each pair as the 
total number of young fledged from all nesting attempts 
in a given year. 

Before analyzing the fecundity data, we visually 
inspected spatial plots and variograms to evaluate the 
degree of spatial autocorrelation in fecundity across our 
study area. A bubble plot of fecundity was created using 
the ggplot2 package in R (Wickham 2009; R Development 
Core Team 2012), and showed no evidence of variation 
in the distribution of fecundity across space (Fig. SI). In 
addition, variograms produced by combining data from 
all years (Fig. S2), as well as those based on each year 
alone, showed little to no increase in the semivariance 



with increasing distance. Variograms were produced 
using the geoR package in R (Diggle and Ribeiro 2007). 
Because these plots showed only weak evidence of spatial 
autocorrelation in fecundity, we proceeded with analyses 
assuming independence between pairs in neighboring 
territories and between pairs breeding in similar locations 
in subsequent years. 

A major goal of our study was to identify the processes 
via which population density and rainfall affected fecun- 
dity. To do so, we used zero-inflated models, which 
explicitly modeled fecundity as a mixture of two probabil- 
ity distributions: one describing the probability of failure 
and one describing the number of offspring produced. 
Zero-inflated models appropriately handle data with a 
surplus of zeros (Lambert 1992) and so are increasingly 
used to analyze patterns of distribution and abundance 
(Fletcher et al. 2005; Martin et al. 2005; Zuur et al. 
2012). Yet despite the potential value of these models for 
demographic analyses of populations in which many indi- 
viduals either do not have the opportunity to breed or 
experience reproductive failure, few studies have used 
these methods to separate the ecological factors that affect 
offspring number from those that affect reproductive fail- 
ure (Quintero et al. 2007; Walker et al. 2009; Smith et al. 
2010). 

In our study population, 52% (95 of 181) of pairs 
fledged no young in a given breeding season, and we used 
zero-inflated Poisson regression to study the sources of 
this demographic variation. While zeros can arise from 
the count-side process, the zero-side process models 
effects that generate a higher proportion of zeros than 
would be expected from a given Poisson distribution. For 
example, a Poisson distribution with a mean of 2.9 (the 
mean number of young fledged by successful pairs; see 
below) should yield data with 5.5% zeros, whereas over 
half of our observations were zeros. The probability of a 
zero was modeled as a binomial distribution with a logit 
link, while the count was modeled as a Poisson distribu- 
tion with a log link. 

We evaluated how rainfall and breeding density (num- 
ber of territorial pairs per ha) affected both the zero-side 
and count-side processes. We calculated breeding season 
precipitation as the total rainfall from November through 
April recorded at a weather station at Middle Ranch, Cat- 
alina Island (33°21'N, 118°26'W). This measure of precip- 
itation was positively correlated with the primary food 
source for the warblers (insect larvae), clutch size, and 
breeding season length, and negatively correlated with 
first clutch initiation date (Sofaer et al. 2013). To account 
for the shared density and rainfall values for pairs in the 
same year, we initially fit a normally distributed random 
effect of year on both the zero and count sides of the 
model. However, the count-side random effect was 
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estimated at zero in the full model, so we refit all models 
with a random effect only on the zero side. We consid- 
ered the following four fixed-effect model structures: 
intercept only, a rainfall effect, a breeding density effect, 
and additive effects of rainfall and breeding density. We 
fit all combinations of these four model structures on 
both the count-side and zero-side processes, for a total of 
16 possible models (Table 1), but did not have sufficient 
years of data to consider interactions between breeding 
density and rainfall. For model selection, we used Ak- 
aike's information criteria adjusted for small sample sizes 
(AICc; (Burnham and Anderson 2002). Fecundity models 
were fit using the NLMIXED procedure in SAS (SAS 
Institute 2008). For all analyses, we report parameter esti- 
mates, standard errors, and confidence intervals based on 
the top model in which each covariate appeared. 

Our zero-inflated models indicated that breeding den- 
sity affected the probability of fledging no young (see 
Results); therefore, to distinguish the effects of population 
density on breeding and/or renesting opportunities from 
its effects on nest predation, we examined how breeding 
density, precipitation, and other factors were correlated 
with daily nest survival probability. We limited this analy- 
sis to nests that were either successful or depredated 
because we were interested in factors affecting nest preda- 
tion, rather than other sources of nest failure (i.e., aban- 
donment, weather, and starvation), and because nest 
predation accounted for 80% of failures once at least one 

Table 1. Model selection results for zero-inflated mixed models of 
fecundity indicated strong support for tlie effects of breeding density 
(bd) and precipitation (precip) on tine probability of fledging zero 
young. We considered all possible additive model structures on both 
the count side and zero side, including intercept-only (.) fixed-effect 
structures. All models contained a normally distributed random effect 
of year on the zero side. 



Count-side 


Zero-side 








-2log 




model 


model 


AICc 


AAlCc 


Weight 


(L) 


k 


bd 


bd + precip 


464.02 


0 


0.33 


451.53 


6 




bd + precip 


464.50 


0.48 


0.26 


454.16 


5 


bd + precip 


bd + precip 


464.58 


0.56 


0.25 


449.93 


7 


precip 


bd + precip 


466.24 


2.22 


0.11 


453.76 


6 


bd 


bd 


471.63 


7.61 


0.01 


461.29 


5 


bd + precip 


bd 


471.65 


7.63 


0.01 


459.17 


6 




bd 


472.06 


8.04 


0.01 


463.84 


4 


bd 


precip 


472.25 


8.23 


0.01 


461.90 


5 


bd + precip 


precip 


472.98 


8.96 


0 


460.50 


6 




precip 


473.34 


9.32 


0 


465.11 


4 


precip 


bd 


473.57 


9.55 


0 


463.23 


5 


bd 




474.84 


10.82 


0 


466.61 


4 


bd + precip 




475.07 


11.05 


0 


464.72 


5 


precip 


precip 


475.15 


11.13 


0 


464.81 


5 






475.54 


11.52 


0 


469.41 


3 


precip 




477.16 


13.14 


0 


468.93 


4 



egg had been laid. We checked nest contents and 
observed adult behavior to assign a nest fate to all nests 
included in our analysis. For each breeding pair in each 
year, we included the earliest successful or depredated 
nest in which at least one egg was laid {n = 142 nests; 
mean per year = 24 nests); this restriction was imposed 
because nests from the same pair in the same year would 
not be independent. Following Rotella et al. (2004, 2007), 
we modeled daily nest survival using a generalized mixed 
model assuming a binomial distribution. Each day a nest 
was known to be active and each interval over which a 
nest failed contributed one survival interval to our analy- 
sis (n = 1983 intervals). All models were fit with a logit 
link function and included a normally distributed random 
effect of year because the nests from each year shared the 
same breeding density and level of precipitation. We built 
models including all additive combinations of the follow- 
ing five fixed effects: precipitation, date, nest height 
(Peluc et al. 2008), breeding density, and nearest neighbor 
distance (Table S2). The number of pairs per ha repre- 
sented plot-level density effects, whereas nearest neighbor 
distance was used to assess the effects of local density 
(these two measures of density were not strongly corre- 
lated; r = —0.28). We used AIQ for model selection. All 
nest survival models were fit using the NLMIXED proce- 
dure in SAS (SAS Institute 2008). Standard errors and 
confidence intervals were estimated with the delta method 
(Powell 2007; Cooch and White 2012). 

Adult survival 

We fit Cormack-JoUy-Seber models (Lebreton et al. 1992) 
to estimate annual apparent survival (</)) and recapture (p) 
probabilities of territorial adults (« = 197) from March to 
the following March. We evaluated how apparent survival 
was influenced by conditions across the annual cycle by 
considering models with effects from both the breeding 
and nonbreeding seasons. Warblers breeding on Catalina 
Island largely winter on the adjacent mainland, so our 
models included covariates describing mainland condi- 
tions during the nonbreeding season. We considered all 
possible additive combinations of the following factors 
(described below) on annual apparent survival probability 
(0): sex, breeding density, November-April precipitation, 
winter population density, and winter precipitation on the 
mainland (Table S3). AH models included effects of sex 
and year on the resighting probability, p, and were fit 
using a logit link function in Program MARK (White and 
Burnham 1999) and ranked based on AICc. 

We calculated an index of winter warbler density based 
on data from the Christmas Bird Count (hereafter CBC; 
National Audubon Society 2010), a 1-day count that 
occurs between mid-December and early January. 
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Volunteers follow designated routes within 24-km radius 
count circles and record all birds observed that day; the 
location of each count circle is consistent between years. 
Groups of participants, known as parties, search for birds, 
and effort is measured as the total number of party hours 
(Dunn et al. 2005). We summarized data from count cir- 
cles in the regular wintering range of O. c. sordida, which 
includes coastal habitats in northern Baja California and 
five southern California counties: Los Angeles, Orange, 
San Diego, Santa Barbara, and Ventura (K. Garrett, pers. 
comm.). W^e included count circles where at least one 
orange-crowned warbler was detected in all years between 
the winters of 2003-04 and 2009-10 (n = 19 count cir- 
cles). CBC data did not identify O. celata to subspecies, 
and counts likely included subspecies other than 
O. c. sordida (Dunn and Garrett 1997). Morphological 
differences between subspecies are relatively slight (Gilbert 
et al. 2010), so we assumed that O. celata subspecies were 
ecologically equivalent during the nonbreeding season 
and averaged the number of orange-crowned warblers 
seen per party hour across all counts to generate an index 
of winter density in each year (range 0.32-0.53 mean 
count per party hour). The use of CBC data necessarily 
entails several simplifying assumptions (reviewed in Dunn 
et al. 2005), and our use of the number of birds per party 
hour as an index of density assumed a linear relationship 
between effort and the resulting count (Link and Sauer 
1999). This latter assumption appeared justified because 
the count data used in our analysis showed no evidence 
of reaching an asymptote with increasing effort. 

We included precipitation in the survival models 
because our previous work in southern California demon- 
strated a strong correlation between rainfall and food 
abundance for songbirds (Morrison and Bolger 2002; 
Sofaer et al. 2013). We used yearly November to AprQ 
precipitation totals for Catalina Island, described above, to 
represent conditions during the breeding season. Winter 
precipitation was calculated by averaging the total Novem- 
ber to February rainfall at Western Regional Climate Cen- 
ter (www.wrcc.dri.edu) weather stations (n = 14) that had 
complete data and were located in the same five California 
counties as the CBC data. Wintering densities on the 
mainland were correlated with mainland precipitation 
during the survival interval (r = —0.75), so we did not 
build models that included both of these covariates. 

Results 

Breeding density and precipitation 

To complement our demographic analyses, we examined 
the relationship between breeding density and precipitation 
to assess how breeding density responded to conditions in 



previous years. Breeding density was not strongly corre- 
lated with total November to April precipitation either in 
the same year (r = —0.13; Fig. S3A) or in the previous year 
(r = 0.26; Fig. S3B). In addition, our exploratory analyses 
were suggestive of population regulation, as breeding den- 
sity in a given year was negatively correlated with density 
in the previous year (r = —0.79; Fig. S4A), perhaps reflect- 
ing the positive relationship between density and mean 
fecundity in the previous year (r = 0.76; Fig. S4B). 

Fecundity and nest survival 

Mean annual fecundity (±1 SE) was 1.40 ± 0.12 young 
fledged per pair per breeding season (n = 181 pairs). 
Excluding 2007, when no birds bred successfully, resulted 
in a mean annual fecundity estimate of 1.55 ± 0.13 
young per pair {n = 163 pairs); the proportion of pairs 
fledging no young in each of these years ranged from 
0.00 to 0.87 (Table SI). Birds that successfijUy fledged at 
least one offspring had a mean annual fecundity of 
2.94 ± 0.11 young per pair (n = 86 pairs). Zero-inflated 
models found weak effects of rainfall and density on the 
number of offspring fledged (Fig. 2A and B), but strong 
effects of these factors on the probability of fledging no 
offspring (Fig. 2C and D). Birds were more likely to be 
unsuccessful in drier years (Fig. 2C; /? = —0.07 ± 0.02; 
95% CI: —0.11, —0.03) and in years with higher breeding 
density (Fig. 2D; /? = 1.64 ± 0.45; 95% CI: 0.53, 2.75). 
Removing precipitation and breeding density from the 
top model increased the AICc value by 7.61 and 8.23, 
respectively, greatly reducing model fit (Table 1). We 
found less support for effects of rainfall and breeding 
density on the count side of the model, that is, on the 
number of offspring fledged (Table 1). The top model 
contained a breeding density effect on the count side, 
but a model without this effect received equivalent sup- 
port (AAICc = 0.48; Table 1). Furthermore, the confi- 
dence interval for the count-side breeding density 
parameter (±1 SE) included zero (/? = -0.14 ± 0.09; 
95% CI: -0.36, 0.07), and support for this effect 
appeared to be influenced by 1 year of data (Fig. 2B). 
Similarly, a model including an effect of precipitation on 
the count side also received equivalent support 
(AAICc = 0.56; Table 1), and here too the confidence 
interval for this parameter included zero (/J = 0.006 
± 0.005; 95% CI: -0.006, 0.018). We therefore con- 
cluded that while support for effects of breeding density 
and rainfall on the probability of reproductive failure was 
strong, we had only weak evidence suggesting these 
variables influenced the number of young fledged from 
successful nests during our study. 

Nest survival showed a strong pattern of density depen- 
dence. Daily survival probability was lower in years with a 



2742 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



H. R. Sofaer et al. 



Partitioning Demographic Variation 



(A) 



4.0 



UJ 

w 
, +1 



.p ^ 



ra cn 



3.5 



3.0 



2.5 



I i 



(B) 



4.0 



3.5 



3.0 



2.5 



I 

A 

i 

T 

A 



Figure 2. The fecundity of pairs that 
successfully fledged young was (A) not 
affected by precipitation and (B) weakly 
affected by breeding density. The proportion 
of pairs fledging no young in each year was 
strongly affected by both (C) precipitation and 
(D) breeding density. In (C), the y-axis 
coordinates are the residuals from a simple 
linear regression of the effects of density on 
the proportion of unsuccessful pairs, while in 
(D), they are the residuals from a simple linear 
regression of the effects of rainfall on the 
proportion of unsuccessful pairs. See text for 
formal zero-inflated analysis of fecundity. 
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higher breeding density (Fig. 3, plotted as daily nest pre- 
dation rate = 1 — daily nest survival rate; = —0.67 
± 0.18; 95% CI: -1.13, -0.20); removing breeding den- 
sity from the top model increased AICc by 7.45 (Table 
S2). Model selection indicated less support for the effects 
of nearest neighbor distance on nest survival (Table S2; 
P = -0.02 ± 0.02; 95% CI: -0.06, 0.02). We found little 
support for an effect of rainfall on nest survival probability 
(Table S2; j3 = -0.01 ± 0.01; 95% CI: -0.03, 0.01). 
Nest survival also declined as the season progressed 
(/? = -0.03 ± 0.01; 95% CI: -0.05, -0.01) and increased 
with nest height, although zero was included near 
the boundary of the 95% confidence interval 
(j3 = 0.20 ± 0.11; 95% CI: -0.09, 0.50). Daily nest sur- 
vival probability (±1 SE) was 0.974 ± 0.004 based on our 
top model and mean breeding density, date, and nest 
height values. To illustrate the magnitude of the effect of 
breeding density on the probability of nest predation, we 
raised the estimated daily nest survival rate to the 12th 
power (the length of the incubation period), and predicted 
90.3% versus 50.7% nest survival over that period at the 
lowest versus the highest observed breeding densities. 

Adult survival 

We found little evidence that precipitation or population 
density during either the breeding or wintering seasons 
affected annual adult survival probability, as the model 



with only a sex effect had equivalent statistical support as 
models containing rainfall and density effects (Table S3). 
Our confidence intervals included zero for the effects of 
breeding density (/? = -0.53 ± 0.44; 95% CI: -1.39, 
0.34), wintering density (^ = -1.79 ± 1.68; 95% 
CI: —5.07, 1.50), November-April precipitation prior to 
the survival interval (precip x—l; P = 0.01 ± 0.01; 95% 
CI: —0.00, 0.03), and November-February precipitation 
during the survival interval (precip x; fl = 0.01 ± 0.01; 
95% CI: —0.01, 0.03). Estimated annual survival probabil- 
ity was higher for males (0.68 ± 0.03; 95% CI: 0.62, 0.73) 
than for females (0.56 ± 0.05; 95% CI: 0.47, 0.66). 

Discussion 

We expected that intraspecific competition for food 
would underlie demographic variation in our study popu- 
lation and would reduce offspring number or lower adult 
survival. Like many insular animal populations (Adler 
and Levins 1994; Blondel 2000), orange-crowned warblers 
on Catalina Island breed at high densities (up to 6.3 
pairs/ha versus up to 1.5 pairs/ha in mainland habitats; 
Gilbert et al. 2010), and male warblers exhibit increased 
levels of testosterone and aggression (Horton et al. 2010; 
Yoon et al. 2012). Nevertheless, breeding density had only 
a weak effect on the number of young produced by suc- 
cessful individuals (Fig. 2B) and was not correlated with 
annual adult survival (Table S3). Instead, breeding density 
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Figure 3. Daily nest predation rate increased witli breeding density. 
Between 2003 and 2009, breeding density varied from 3.4 to 6.3 
territories/tia. 

affected the probability of reproductive failure (Fig. 2D) 
via density-dependent nest predation (Fig. 3). By parti- 
tioning the factors affecting reproductive failure from 
those affecting the fecundity of successful individuals, we 
were able to improve our mechanistic understanding of 
demographic variation and separate the effects of density- 
dependent nest predation, intraspecific competition, and 
precipitation. 

Because we found that both breeding density and rain- 
fall affected the probability of fledging zero young 
(Fig. 2C and D), we evaluated the influence of these fac- 
tors on three major processes that determine whether a 
territorial individual wUl reproduce successfully in a given 
year. First, a bird must breed, rather than skip breeding 
for a year. Second, nest survival depends on attributes of 
the parent, nest site, and predator community (reviewed 
in Cresswell 1997). Third, following nest failure, birds 
either renest or suspend breeding and enter the nonbreed- 
ing period, a decision that repeats with each failed repro- 
ductive attempt. The probability of renesting depends on 
the timing of failure and may be affected by food avail- 
ability or the energetic reserves of the breeding adult 
(Nagy and Holmes 2005; Arnold et al. 2010; Caldwell 
et al. 2013). Below, we discuss how the effects of rainfall 
and breeding density may have been mediated via 
these three sources of reproductive failure in O. celata on 
Catalina Island. 

Rainfall likely influenced the probability of reproduc- 
tive failure by altering birds' propensity to initiate breed- 
ing and to renest after nest failure. Indeed, the 
importance of rainfall was observed during a severe 
drought, which led the majority of territorial pairs to skip 
breeding in 2007; not a single pair successfully fledged 



young in that year (Fig. 2C at 9 cm of rainfall; Langin 
et al. 2009). Increased rainfall was also associated with 
longer breeding seasons (Sofaer et al. 2013), suggesting 
that the relationship between rainfall and food abundance 
led birds to increase their probability of renesting as pre- 
cipitation increased. Rainfall was not correlated with daily 
nest survival probability (Table S2). We expected that 
rainfall would have had a stronger effect on the number 
of young fledged (Fig. 2A) because birds in this popula- 
tion lay larger clutches in wetter years (Sofaer et al. 
2013). However, the relatively small magnitude of change 
in mean clutch size between years (range: 2.8-3.5 eggs, 
including all nests throughout the breeding season) could 
explain our failure to find the expected effects on off- 
spring number. 

Nest predation appeared to be the primary mechanism 
underlying the density-dependent increase in reproductive 
failure (Figs. 2D, 3). Breeding density had no apparent 
effects on the probability of breeding, as all territorial 
pairs attempted to breed in the year with the highest 
observed conspecific density. Breeding density also did 
not influence breeding season length; a post hoc analysis 
of breeding season length found that the addition of 
breeding density to a simple linear model including pre- 
cipitation increased the AIC value by 1.74, and the confi- 
dence interval on the density effect included zero. 

Density-dependent nest predation could arise from 
either functional or numeric responses of nest predators 
to increasing prey density (HoUing 1959) or from a form 
of site dependence in which areas with lower nest preda- 
tion risk are preferentially occupied (Rodenhouse et al. 
2003; Emmering and Schmidt 2011). This latter mecha- 
nism could arise if warblers can accurately assess spatial 
variation in nest predation risk, but testing this hypothe- 
sis is difficult because the probability of nest predation 
risk may shift spatially between years so particular territo- 
ries may not consistently be favored. An alternative expla- 
nation of our results, that the relationship between 
warbler density and nest predation rates simply reflected 
correlated numerical responses of warbler and nest preda- 
tor populations to high productivity in wet years, is unli- 
kely because warbler densities were not correlated with 
precipitation in the previous year (Fig. S3B). Nevertheless, 
tracking nest predator abundance would allow for an 
analysis of the relationship between predator density 
and nest predation rates, which remains understudied 
(Abrams and Ginzburg 2000; Oro et al. 2006; Schmidt 
et al. 2008). 

There is a need to understand how reptilian and mam- 
malian predators locate nests and allocate foraging effort 
(Weatherhead and Bloun-Demers 2004), particularly 
given the possibility that one or more of the island's 
predators may not simply find nests incidentally. 



2744 



© 2014 The Autliors. Ecology and Evolution publisiied by John Wiley & Sons Ltd. 



H. R. Sofaer et al. 



Partitioning Demographic Variation 



Documented or likely nest predators in this system 
include gopher snakes (Pituophis catenifer), California 
kingsnakes (Lampropeltis calif orniae), island foxes {Urocy- 
on littoralis), feral cats {Felis catus), deer mice {Peromyscus 
maniculatus), and black and Norwegian rats Rattus rattus, 
R. norvegicus; Catalina Island lacks avian nest predators 
such as jays (Peluc et al. 2008). Currently, there is no 
consensus about how differences in the type, diversity, 
and abundance of predator and prey species influence the 
likelihood of nest predation being density dependent, 
which has been found in several systems (Andersson and 
Wiklund 1978; Martin 1996; Lariviere and Messier 1998; 
Schmidt and Whelan 1999; Gunnarsson and Elmberg 
2008) but not in others (Zimmerman 1984; O'Reilly and 
Hannon 1989; Reitsma 1992; Ackerman et al. 2004; Sillett 
and Holmes 2005). However, our study is one of the few 
that has both documented density-dependent nest preda- 
tion in a natural population and also assessed whether 
nest predation underlies variation in fecundity, giving it 
the potential to regulate the focal population (Arcese 
et al. 1992; Tapper et al. 1996). 

Density-dependent nest predation also has implications 
for the evolution and expression of nest-site selection. 
When a nest predator's functional response includes spe- 
cializing on a set of nest-site characteristics, individuals 
may benefit by nesting in less-used sites, a process that 
can affect community composition and the evolution of 
nest-site selection (Martin 1988, 1996), although not all 
studies have found support for predator specialization 
(Reitsma and Whelan 2000; Rangen et al. 2001). Orange- 
crowned warblers on Catalina Island (O. c. sordida) exhi- 
bit an unusual amount of variation in nest height and 
nest-site location relative to most other Oreothlypis, 
including mainland populations of O. celata, which 
almost always nest on the ground. Warbler nest sites on 
Catalina Island range from ground nests that experience 
relatively high nest predation, to safer off-ground nests in 
shrubs and tree crowns (Peluc et al. 2008). Our results 
therefore raise the possibility that density-dependent nest 
predation could favor the diversification of nesting sites 
within a population, rather than solely between species. 
Additional research is needed to evaluate nest predators' 
search methods and functional responses, whether high 
warbler breeding densities on Catalina Island could have 
favored specialization by predators, and whether warblers 
gain any selective advantage from choosing atypical nest- 
site locations; support for this latter hypothesis would 
imply that nest predation is frequency dependent. 

Conclusions 

Density-dependent nest predation appeared to be an 
important regulatory mechanism for O. celata on Catalina 



Island (Fig. 3), even though our study population bred at 
high densities and exhibited traits associated with strong 
intraspecific competition (Horton et al. 2010; Yoon et al. 
2012). However, our research focused only on territory 
holders, and intraspecific competition for food and terri- 
tories is likely to influence postfledging survival and 
recruitment in this population. Our ability to interpret 
patterns of density-dependent fecundity was facilitated by 
the use of zero-inflated models, which allowed us to sepa- 
rate the ecological factors affecting the count-side process 
(i.e., the number of young fledged) from those affecting 
the zero-side process (i.e., the probability of fledging 
no young). We suggest that the broader use of zero- 
inflated models will improve ecological inference about 
populations in which many individuals fail to successfully 
reproduce. 

Acknowledgments 

This work was supported by The Nature Conservancy, 
the Smithsonian Institution, and Colorado State Univer- 
sity. HRS was supported by Colorado State University, 
NSF-IGERT Grant DGE-#0221595 (administered by the 
PRIMES program at Colorado State University), and 
National Science Foundation DEB-0846175 to CKG. 
Research was conducted under Colorado State University 
lACUC permit 06-043A-01 and Smithsonian Institution 
permits # 06-03 and 07-24, and US federal banding 
permit 22665 issued to the Smithsonian Migratory Bird 
Center. The Catalina Island Conservancy provided logisti- 
cal support. We thank the many field biologists who 
helped with data collection, particularly J. Coumoutso, H. 
Montag, S. Peluc, and J. Yoon. V. Bakker, P. Chapman, 
and J. Rotella provided statistical advice. This paper 
benefited from comments from anonymous reviewers, C. 
Handelsman, B. Noon, E. Ruell, C. Webb, and L. Zanette. 

Data Accessibility 

Data are available from the Dryad Digital Repository 
doi:10.5061/dryad.m4fs7. 

Conflict of Interest 

None declared. 
References 

Abrams, P. A., and L. R. Ginzburg. 2000. The nature of 

predation: prey dependent, ratio dependent or neither? 

Trends Ecol. Evol. 15:337-341. 
Ackerman, I. T., A. L. Blackmer, and ]. M. Eadie. 2004. Is 

predation on waterfowl nests density dependent? Tests at 

three spatial scales. Oikos 107:128-140. 



2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



2745 



Partitioning Demographic Variation 



H. R. Sofaer ef al. 



Adler, G. H., and R. Levins. 1994. The island syndrome in 

rodent populations. Q. Rev. Biol. 69:473-490. 
Andersson, M., and C. G. Wiklund. 1978. Clumping versus 

spacing out - experiments on nest predation in fieldfares 

{Turdus pilaris). Anim. Behav. 26:1207-1212. 
Arcese, P., and J. N. M. Smitti. 1988. Effects of population 

density and supplemental food on reproduction in song 

sparrows. J. Anim. Ecol. 57:119-136. 
Arcese, P., J. N. M. Smith, W. M. Hochachka, C. M. Rogers, 

and D. Ludwig. 1992. Stability, regulation, and the 

determination of abundance in an insular song sparrow 

population. Ecology 73:805-822. 
Arnold, T. W., J. H. Devries, and D. W. Howerter. 2010. 

Factors that affect renesting in mallards (Anas 

platyrhynchos). Auk 127:212-221. 
Blondel, J. 2000. Evolution and ecology of birds on islands: 

trends and prospects. Vie Et Milieu-Life Environ. 50:205-220. 
Brouwer, L., J. M. Tinbergen, C. Both, R. Bristol, D. S. 

Richardson, and J. Komdeur. 2009. Experimental evidence 

for density-dependent reproduction in a cooperatively 

breeding passerine. Ecology 90:729-741. 
Burnham, K. P., and D. R. Anderson. 2002. Model selection 

and multimodel inference: a practical information-theoretic 

approach. Springer, New York. 
Caldwell, L., V. J. Bakker, T. S. Sillett, M. A. Desrosiers, S. A. 

Morrison, and L. M. Angeloni. 2013. The influence of 

habitat and nest observation on the reproductive success of 

the island scrub-jay. Condor 115:603-613. 
Glutton-Brock, T. H. 1988. Reproductive success. The 

University of Chicago Press, Chicago. 
Glutton-Brock, T. H., M. Major, S. D. Albon, and F. E. 

Guinness. 1987. Early development and population 

dynamics in red deer. I. Density-dependent effects on 

juvenile survival. J. Anim. Ecol. 56:53-67. 
Cooch, E. G., and G. C. White. 2012. Program MARK: a 

gentle introduction [Online]. 11th edn. http://www.phidot. 

org/software/mark/docs/book/ (accessed 12 April 2012). 
Coulson, T., E. A. Catchpole, S. D. Albon, B. J. T. Morgan, 

J. M. Pemberton, T. H. Glutton-Brock, et al. 2001. Age, sex, 

density, winter weather, and population crashes in Soay 

sheep. Science 292:1528-1531. 
Covas, R. 2012. Evolution of reproductive life histories in 

island birds worldwide. Proc. Biol. Sci. 279:1531-1537. 
Cresswell, W. 1997. Nest predation: the relative effects of nest 

characteristics, clutch size and parental behaviour. Anim. 

Behav. 53:93-103. 
Dempster, J. P. 1971. Population ecology of cinnabar moth, 

Tyria jacobaeae L. (Lepidoptera, Arctiidae). Oecologia 

7:26-67. 

Development Core Team, R. 2012. R: A language and 
environment for statistical computing. R Foundation for 
Statistical Computing, Vienna, Austria. ISBN 3-900051-07-0. 

Diggle, P. J., and P. J. Ribeiro. 2007. Model-based geostatistics. 
Springer, New York, NY. 



Dunn, J., and K. Garrett. 1997. A field guide to warblers of 
North America. Houghton Mifflin Harcourt, Boston. 

Dunn, E. H., C. M. Francis, P. J. Blancher, S. R. Drennan, M. 
A. Howe, D. Lepage, et al. 2005. Enhancing the scientific 
value of The Christmas Bird Count. Auk 122:338-346. 

Eggers, S., M. Griesser, M. Nystrand, and J. Ekman. 2006. 
Predation risk induces changes in nest-site selection and 
clutch size in the Siberian jay. Proc. Biol. Sci. 273:701-706. 

Emmering, Q. C., and K. A. Schmidt. 2011. Nesting songbirds 
assess spatial heterogeneity of predatory chipmunks by 
eavesdropping on their vocalizations. J. Anim. Ecol. 
80:1305-1312. 

Fletcher, D., D. MacKenzie, and E. Villouta. 2005. Modelling 
skewed data with many zeros: a simple approach combining 
ordinary and logistic regression. Environ. Ecol. Stat. 
12:45-54. 

Forchhammer, M. C, N. C. Stenseth, E. Post, and R. 

Langvatn. 1998. Population dynamics of Norwegian red 

deer: density-dependence and climatic variation. Proc. Biol. 

Sci. 265:341-350. 
Forrester, G. E. 1995. Strong density-dependent survival and 

recruitment regulate the abundance of a coral reef fish. 

Oecologia 103:275-282. 
Frederiksen, M., and T. Bregnballe. 2000. Evidence for 

density-dependent sur'vival in adult cormorants from a 

combined analysis of recoveries and resightings. J. Anim. 

Ecol. 69:737-752. 
GaiUard, J. M., M. Festa-Bianchet, and N. G. Yoccoz. 1998. 

Population dynamics of large herbivores: variable 

recruitment with constant adult survival. Trends Ecol. Evol. 

13:58-63. 

Gilbert, W. M., M. K. Sogge, and C. Ill Van Riper. 2010 
Orange-crowned Warbler [Oreothlypis celata). A. Poole, ed. 
The birds of North America online. Cornell Lab of 
Ornithology, Ithaca. Retrieved from the Birds of North 
America Online: http://bna.birds.corneU.edu/bna/species/101 
doi:10.2173/bna.l01 (accessed 4 November 2013). 

Godfray, H. C. J., L. Partridge, and P. H. Har-vey. 1991. Clutch 
size. Annu. Rev. Ecol. Syst. 22:409-429. 

Grant, P. R., B. R. Grant, L. F. Keller, and K. Petren. 2000. 
Effects of El Nino events on Darwin's finch productivity. 
Ecology 81:2442-2457. 

Gunnarsson, G., and J. Elmberg. 2008. Density-dependent nest 
predation - an experiment with simulated Mallard nests in 
contrasting landscapes. Ibis 150:259-269. 

Harms, K. E., S. J. Wright, O. Calderon, A. Hernandez, and 
E. A. Herre. 2000. Pervasive density-dependent recruitment 
enhances seedling diversity in a tropical forest. Nature 
404:493-495. 

Hrxon, M. A., and M. H. Carr. 1997. Synergistic predation, 
density dependence, and population regulation in marine 
fish. Science 277:946-949. 

Hochachka, W. M., and A. A. Dhondt. 2000. 

Density-dependent dechne of host abundance resulting from 



2746 



© 2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



H. R. Sofaer et al. 



Partitioning Demographic Variation 



a new infectious disease. Proc. Natl Acad. Sci. USA 97: 
5303-5306. 

Holbrook, S. J., and R. J. Sclimitt. 2002. Competition for 
shelter space causes density-dependent predation mortality 
in damselfishes. Ecology 83:2855-2868. 

HoUing, C. 1959. The components of predation as revealed by 
a study of small-mammal predation of the European pine 
sawfly. Can. Entomol. 91:293-320. 

Horton, B. M., J. Yoon, C. K. Ghalambor, I. T. Moore, and 
T. S. SiUett. 2010. Seasonal and population variation in male 
testosterone levels in breeding orange-crowned warblers 
(Vermivora celata). Gen. Comp. Endocrinol. 168:333-339. 

Jenouvrier, S., C. Barbraud, and H. Weimerskirch. 2003. 
Effects of climate variability on the temporal population 
dynamics of southern fulmars. J. Anim. Ecol. 72:576-587. 

Lack, D. 1947. The significance of clutch size. Ibis 89:302-352. 

Lambert, D. 1992. Zero-inflated Poisson regression, with an 
application to defects in manufacturing. Technometrics 
34:1-14. 

Langin, K. M., T. S. SiUett, J. Yoon, H. R. Sofaer, and C. K. 
Ghalambor. 2009. Reproductive consequences of an extreme 
drought for songbirds on Santa Catalina and Santa Cruz 
Islands. Pp. 293-300 in C. C. Damiani and D. K. Garcelon, 
eds. Proceedings of the seventh California Islands 
symposium. Institute for Wildlife Studies, Areata. 

Lariviere, S., and F. Messier. 1998. Effect of density and nearest 
neighbours on simulated waterfowl nests: can predators 
recognize high-density nesting patches? Oikos 83:12-20. 

Lebreton, J. D., K. P. Burnham, J. Clobert, and D. R. 
Anderson. 1992. Modeling survival and testing biological 
hypotheses using marked animals - a unified approach with 
case-studies. Ecol. Monogr. 62:67-118. 

Link, W. A., and J. R. Sauer. 1999. Controlling for varying 
effort in count surveys - An analysis of Christmas Bird 
Count data. J. Agric. Biol. Environ. Stat. 4:116-125. 

MacArthur, R. H., J. M. Diamond, and J. R. Karr. 1972. 

Density compensation in island faunas. Ecology 53:330-342. 

Martin, T. E. 1987. Food as a limit on breeding birds - a 
life-history perspective. Annu. Rev. Ecol. Syst. 18:453-487. 

Martin, T. E. 1988. On the advantage of being different - nest 
predation and the coexistence of bird species. Proc. Natl 
Acad. Sci. USA 85:2196-2199. 

Martin, T. E. 1996. Fitness costs of resource overlap among 
coexisting bird species. Nature 380:338-340. 

Martin, T. E., P. R. Martin, C. R. Olson, B. J. Heidinger, and 
J. J. Fontaine. 2000. Parental care and clutch sizes in North 
and South American birds. Science 287:1482-1485. 

Martin, T. G., B. A. Wintie, J. R. Rhodes, P. M. Kuhnert, S. A. 
Field, S. J. Low-Choy, et al. 2005. Zero tolerance ecology: 
improving ecological inference by modelling the source of 
zero observations. Ecol. Lett. 8:1235-1246. 

Morrison, S. A., and D. T. Bolger. 2002. Variation in a 
sparrow's reproductive success with rainfall: food and 
predator-mediated processes. Oecologia 133:315—324. 



Murdoch, W. W. 1994. Population regulation in theory and 

practice. Ecology 75:271-287. 
Nagy, L. R., and R. T. Holmes. 2005. Food limits annual 

fecundity of a migratory songbird: an experimental study. 

Ecology 86:675-681. 
National Audubon Society. 2010. The Christmas Bird Count 

Historical Results [Online]. Available at http://www. 

christmasbirdcount.org. (accessed 1 February 2012). 
Nevoux, M., O. Gimenez, D. Arlt, M. NicoU, C. Jones, and 

K. Norris. 2011. Population regulation of territorial species: 

both site dependence and interference mechanisms matter. 

Proc. Biol. Sci. 278:2173-2181. 
Newton, I. 1998. Population limitation in birds. Academic 

Press, London. 

O'Reilly, P., and S. J. Hannon. 1989. Predation of simulated 
wiUow ptarmigan nests - the influence of density and cover 
on spatial and temporal patterns of predation. Can. J. Zool. 
67:1263-1267. 

Oro, D., A. Martinez-Abrain, M. Paracuellos, J. C. Nevado, 

and M. Genovart. 2006. Influence of density dependence on 

predator-prey seabird interactions at large spatio-temporal 

scales. Proc. Biol. Sci. 273:379-383. 
Ostfeld, R. S., and C. D. Canham. 1995. Density-dependent 

processes in meadow voles - an experimental approach. 

Ecology 76:521-532. 
Peluc, S. I., T. S. SiUett, J. T. Rotenberry, and C. K. Ghalambor. 

2008. Adaptive phenotypic plasticity in an island songbird 

exposed to a novel predation risk. Behav. Ecol. 19:830-835. 
Powell, L. A. 2007. Approximating variance of demographic 

parameters using the delta method: a reference for avian 

biologists. Condor 109:949-954. 
Quintero, H. E., A. Abebe, and D. A. Davis. 2007. 

Zero-inflated discrete statistical models for fecundity data 

analysis in channel catfish, Ictalurus punctatus. J. World 

Aquaculture Soc. 38:175-187. 
Rangen, S. A., R. G. Clark, and K. A. Hobson. 2001. Predator 

responses to similarity and dispersion of artificial nest sites: 

implications for the structure of Boreal forest songbird 

communities. Auk 118:105-115. 
Reitsma, L. 1992. Is nest predation density dependent - a test 

using artificial nests. Can. J. Zool. 70:2498-2500. 
Reitsma, L. R., and C. J. Whelan. 2000. Does vertical 

partitioning of nest sites decrease nest predation? Auk 

117:409-415. 

Ricklefs, R. E. 1969. An analysis of nesting mortality in birds. 

Smithson. Contrib. Zool. 9:1-48. 
Rodenhouse, N. L., T. S. SiUett, P. J. Doran, and R. T. 

Holmes. 2003. Multiple density-dependence mechanisms 

regulate a migratory bird population during the breeding 

season. Proc. Biol. Sci. 270:2105-2110. 
RoteUa, J., S. Dinsmore, and T. Shaffer. 2004. ModeUng 

nest-survival data: a comparison of recently developed 

methods that can be implemented in MARK and SAS. 

Anim. Biodivers. Conserv. 27:187-205. 



2014 The Authors. Ecology and Evolution published by John Wiley & Sons Ltd. 



2747 



Partitioning Demographic Variation 

Rotella, J., M. Taper, S. Stephens, and M. Lindberg. 2007. 
Extending methods for modehng heterogeneity in 
nest-survival data using generalized mixed models. Stud. 
Avian Biol. 34:34-44. 

SAS Institute. 2008. SAS version 9.2. SAS Institute, Gary, NC. 

Schmidt, K. A., and C. J. Whelan. 1999. Nest predation on 
woodland songbirds: when is nest predation density 
dependent? Oikos 87:65-74. 

Schmidt, K. A., S. A. Rush, and R. S. Ostfeld. 2008. Wood 
thrush nest success and post-fledging survival across a 
temporal pulse of small mammal abundance in an oak 
forest. J. Anim. Ecol. 77:830-837. 

SiUett, T. S., and R. T. Holmes. 2005. Long-term demographic 
trends, limiting factors, and the strength of density 
dependence in a breeding population of a migratory 
songbird. Pp. 426-436 in R. S. Greenberg and P. P. Marra, 
eds. Birds of two worlds: the ecology and evolution of 
temperate-tropical migration. Johns Hopkins University 
Press, Baltimore. 

SiUett, T. S., N. L. Rodenhouse, and R. T. Holmes. 2004. 
Experimentally reducing neighbor density affects 
reproduction and behavior of a migratory songbird. Ecology 
85:2467-2477. 

Skutch, A. F. 1949. Do tropical birds rear as many young as 
they can nourish. Ibis 91:430-458. 

Smith, S., T. Mang, J. G. De Bellocq, H. Schaschl, 

C. Zeitlhofer, K. Hacklander, et al. 2010. Homozygosity at a 
class II MHC locus depresses female reproductive ability in 
European brown hares. MoL Ecol. 19:4131^143. 

Sofaer, H. R., T. S. SiUett, S. I. Peluc, S. A. Morrison, and 
C. K. Ghalambor. 2013. Differential effects of food 
availability and nest predation risk on avian reproductive 
strategies. Behav. Ecol. 24:698-707. 

Southern, H. N. 1970. The natural control of a population of 
Tawny owls {Strix aluco). J. Zool. 162:197-285. 

Tapper, S. C., G. R. Potts, and M. H. Brockless. 1996. The 
effect of an experimental reduction in predation pressure on 
the breeding success and population density of grey 
partridges Perdix perdix. J. Appl. Ecol. 33:965-978. 

Walde, S. J., and W. W. Murdoch. 1988. Spatial density 

dependence in parasitoids. Annu. Rev. Entomol. 33:441^66. 

Walker, M., A. HaU, R. M. Anderson, and M. G. Basanez. 
2009. Density-dependent effects on the weight of female 
Ascaris lumbricoides infections of humans and its impact on 
patterns of egg production. Parasit. Vectors 2:11. 

Weatherhead, P. J., and G. Bloun-Demers. 2004. 

Understanding avian nest predation: why ornithologists 
should study snakes. J. Avian Biol. 35:185-190. 



H. R. Sofaer ef al. 

White, G. C., and K. P. Burnham. 1999. Program MARK: 

survival estimation from populations of marked animals. 

Bird Stud. 46:120-139. 
Wickham, H. 2009. ggplot2: elegant graphics for data analysis. 

Springer, New York. 
Yeaton, R. I. 1974. Ecological analysis of chaparral and pine 

forest bird communities on Santa Cruz Island and mainland 

California. Ecology 55:959-973. 
Yoon, J., T. S. SiUett, S. A. Morrison, and C. K. Ghalambor. 

2012. Breeding density, not life history, predicts 

interpopulation differences in territorial aggression in a 

passerine bird. Anim. Behav. 84:515-521. 
Zanette, L. Y., A. F. White, M. C. AUen, and M. Clinchy. 2011. 

Perceived predation risk reduces the number of offspring 

songbirds produce per year. Science 334:1398-1401. 
Zimmerman, J. L. 1984. Nest predation and its relationship 

to habitat and nest density in dickcissels. Condor 86:68- 

72. 

Zuur, A. F., A. A. Saveliev, and E. N. leno. 2012. Zero inflated 
models and generalized linear mixed models with R. 
Highland Statistics Ltd., Newburgh. 

Supporting Information 

Additional Supporting Information may be found in the 
online version of this article: 

Table SI. Summary statistics describing mean annual var- 
iation in fecundity (total offspring fledged per pair). 
Table 82. Our generalized mixed models showed strong 
support for density-dependent effects on daily nest sur- 
vival; at higher breeding densities daily nest survival was 
lower. 

Table S3. Model selection results for the top ten Cor- 
mack-JoUy-Seber mark-recapture models. 
Figure SI. A bubble plot of the annual fecundity of war- 
bler pairs in our study area showed no consistent pattern 
in the distribution of fecundity across space. 
Figure S2. A variogram of fecundity showing weak evi- 
dence of spatial autocorrelation within our study area. 
Figure S3. Breeding density in a given year was not corre- 
lated with total November-April precipitation (A) in the 
same year or (B) in the previous year. 
Figure S4. Breeding density was (A) negatively correlated 
with breeding density in the previous year, and (B) posi- 
tively correlated with mean fecundity in the previous year. 
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